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The cross-over region of the quark-gluon plasma (QGP) created in heavy ion colhsions is influ- 
enced by the nearby deconfinement, chiral and baryon condensation phase transitions. A char- 
acteristic signature of the deconfinement transition in this region can be inferred using the flux 
tube model, which is dual to the Polyakov loop description and which offers a visual picture of 
what happens during the transition. The three-point (anti)vertices of a flux tube network lead to 
formation of (anti)baryons upon hadronisation. Since there is no fundamental interaction associ- 
ated with the baryon number, correlations in the baryon number distribution at the last scattering 
surface directly reflect the preceding pattern of the flux tube vertices in the QGP. An alternating 
pattern of vertices and antivertices would lead to an oscillatory signal in the two-point baryon 
number correlations, under the experimental conditions prevalent in heavy ion collisions at RHIC 
and LHC. The strength of the oscillations is a measure of the flexibility of the QGP. 
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The study of collective non-peiturbative phenomena in QCD at high temperatures and/or at 
large chemical potentials is a highly active area of research: experimentally thi^ough heavy ion col- 
lisions [|]], theoretically through phenomenological models [|2|], and computationally through lattice 
QCD simulations Consider QCD with N = 3 colours and Nf degenerate quark flavours of mass 
m, at temperature T and quark chemical potential jU. The phase structure of this theory is depicted 
in Fig|l](a) in a schematic manner. Three corners of the phase structure are well-understood because 
of the exact symmetries present there: (a) the finite temperature deconfinement phase transition at 
m = oo, governed by the global Z3 centre symmetry of the Polyakov loop, (b) the finite temperature 
chiral phase transition atm = = jj., governed by the flavour SU {Nf)L ® SU {Nf)R symmetry, and 
(c) the baryon condensation phase transition at m = = T, when jj. crosses the constituent quark 
mass. These end-points are first-order phase transitions, and phase transition surfaces extend in- 
wards from them. Phenomenological and numerical studies show that the surfaces end in critical 
lines, and there is no phase transition for the physical values of the quark masses as T is varied (un- 
less n is sufficiently large). Nevertheless, the three nearby transitions cause various QCD properties 
to change rapidly in the cross-over region. To enhance our understanding of QCD, it is important 
to construct observables that can be extracted from the experimental data and that can highlight the 
dynamical features of deconfinement, chiral symmetry restoration and baryon condensation. 

In this article, I describe how the deconfinement phase transition leaves a signature in two- 
point baryon number correlations in the angular distribution of the hadrons produced in heavy ion 
collisions. The collisions at RHIC and LHC create a high energy state of QCD matter, which 
upon cooling produces 0(1000) hadrons. These hadrons are dominated by pions, but also include 
0(100) baryons and antibaryons. The particle numbers are lai^ge enough for one to analyse their 
distributions in position, momenta and energy (instead of only looking at individual particle prop- 
erties). For the statistical distributions to be meaningful, they have to arise from some reasonably 




(a) (b) 

Figure 1: (a) Schematic description of the QCD phase structure in the m—T — ji space. The first-order 
transition surfaces are shown shaded, and the critical lines are shown dotted. The value of m corresponding 
to the real world QCD, and the QGP region accessible in heavy-ion collisions, are pointed out. 
(b) Schematic representation of the anticipated baryon number two-point correlation functions, (r) (thick 
line) and gy{r) (thin line). The former is similar to that for objects with hard-core repulsion. The latter is for 
a percolating flux tube network where vertices and antivertices alternate. 
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equilibrated state. Experimental data suggest that such an approximately thermalised quark-gluon 
plasma (QGP) is created, with energy density e ~ IGeV/fm^ and temperature T^x ~ 175MeV, in 
the fireball of central collisions and produces moderate pt hadrons. An important limitation of the 
experiments is that they detect only charged hadrons moving in sufficiently transverse directions. 

The experimentally observed multiplicities and distributions of various particles have been fit- 
ted to hadron resonance gas models, resulting in estimates of the "chemical freeze-out" temperature 
(where inelastic scattering stops) and the "kinetic freeze-out" temperature (where elastic scattering 
stops) iQ]. The stage is now set to go beyond the average properties and investigate the two-point 
correlations in the distributions. A useful analogy is the study of the cosmic microwave background 
radiation (CMBR), where first the black body spectrum was determined, and then temperature fluc- 
tuations were detected at the level of \T ~ 10^^ T with characteristic angular correlations [Q]. The 
"last scattering surface" in the evolution of the CMBR is analogous to the "kinetic freeze-out" 
stage of the QCD fireball. The patterns that can be observed there would be an invaluable help in 
understanding the QGP dynamics that preceded it. 

In what follows, I first present a flux tube model providing a physical picture of the deconfine- 
ment process ||^, and then point out that under suitable conditions the scenario predicts a specific 
two-point baryon number coiTclation signal [^, illustrated in Fig.[I](b). The challenge for theorists 
is to quantify this signal as much as possible, and the challenge for experimentalists is to detect it. 

1. Flux Tube Model Phenomenology 

The flux tube model of QCD is motivated by the dual superconductor description of linear- 
colour confinement ||8|], where condensation of colour magnetic charges restricts colour electric 
fields to vortex-like configurations. The model is phenomenologically quite successful, and is 
manifest in lattice QCD results as the area law for large Wilson loops. A characteristic property of 
the flux tube is its energy per unit length, i.e. the string tension a. Other than that, the flux tube has 
a finite width w and a persistence length a, both of order Ag^^. 

The flux tubes have to obey the constraint of Gauss's law. So they terminate only on quarks, 
and interact only at 3-point vertices. These two features represent the invariant tensors dab and Sabc 
used to describe the meson and the baryon wavef unctions, and are illustrated by the lattice QCD 
simulation results displayed in Fig.^ Other multi-quark hadron states are phenomenologically not 
prominent, and so all other interactions among the flux tubes are ignored in the model. 




Figure 2: Gluon flux tube configurations between static (anti)quark sources determined by lattice QCD 
simulations [^. A meson is on the left and a baryon is on the right. 
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The finite temperature behavior of the model is governed by competition between the energy 
and the entropy of the flux tube configurations. At low temperatures, energy wins, keeping down 
the total length of the flux tubes. At high temperatures, entropy dominates, producing elaborate 
structures of the flux tubes all over the space. As the temperature is increased the flux tubes oscillate 
more, and also produce more vertices. Some of the possible configurations are shown in Fig.^a)- 
(d). First consider the pure gauge theory, i.e. m = oo. In absence of vertices, there is a second-order 
deconfinement phase transition, where the flux tube length diverges and the quark-antiquark pair 
loses information about each-other's position. In presence of vertices, the flux tubes can percolate 
the space in a network before the correlation length diverges. That also allows the quark-antiquark 
pair (hooked on to the network) to lose information about each-other's position, but produces a 
first-order deconfinement phase transition. 

When finite mass quarks are included in the model, they break the flux tubes by quark- 
antiquark pair production from the vacuum. With the flux tube network breaking up, the decon- 
finement phase transition weakens as the quark mass is lowered from m = oo and/or as the chemical 
potential is increased. Numerical estimates show that for Nf = 3 QCD at small chemical potentials, 
the first-order deconfinement phase transition ends in a critical line around m = 1 .5GeV, as depicted 
in Fig|l](a). Although the cross-over region does not have any sharp signal of the deconfinement 
phase transition, the percolating flux tube scenario still suggests a detectable observable there. 

It is straightforwai^d to formulate the model on a lattice with spacing a. The flux tubes live on 
the links, while the quarks and vertices live on the sites. The variables take values 0, ±1, depending 
on the direction of the flux, as shown in Fig.^e). The total energy of a flux tube configuration is 

E = aaj^\ni^^ \ +m'^\pij\ +v'^\qi\ , (1.1) 

i.H i.f i 

where v denotes the energy cost of a vertex, and / sums over the 2Nf spin and flavour quark degrees 
of freedom. The constraint of Gauss's law at every site is 




Figure 3: (a)-(d) depict possible flux tube configurations connecting a static quark-antiquark pair. Increasing 
temperature takes (a) to (b) and (c) to (d), while creation of baryonic vertices takes (a) to (c) and (b) to (d). 
(e) lists the link and the site variables used for formulating the flux tube model on a lattice. 
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The grand canonical partition function for the system, with baryon number B, is 



Z[T,^\= £ exp 



E-3IJ.B 



B 



1 



(1.3) 



It is fully factorised by expressing the constraint at every site as 5a,,o = J^ni'^^i/'^^)^'"'^' ■ The 
sum over the variables n,^, pjj, qi can then be carried out explicitly, resulting in 



n(l + 2e"'"/^cos(0, + / 



n(l+2e"'/^cos(30O) • (1-4) 



This is in the universality class of the 3-dim XY spin model, in the presence of an ordinary magnetic 
field as well as a Z3 symmetric magnetic field. In this form, Z[r,/x] possesses a generalised 
symmetry, with ^ corresponding to 0; — —0,- and ^ corresponding to complex conjugation. As a 
consequence, even though Eq.( |1.4| ) involves complex weights, the equivalent form in Eq.(1.3) can 



be numerically simulated without any fermion sign problem at finite chemical potential [ IC , 1 1 1. 

Physical interpretation of the site variables 0, is uncovered by insertion of static sources in the 
system. A static quark at site j modifies the Gauss's law constraint there as — )■ Saj.-i, and 
its free energy is given by exp{—Fq/T) = (exp(— In the gauge field theory language, this 
quantity is the expectation value (Pj) of the Polyakov loop at site j. So we infer that dj represents 
the phase of the Polyakov loop Pj, and the flux tube description of deconfinement is dual to the 
familiar Polyakov loop description of deconfinement. The benefit of the flux tube description is the 
visual representation it provides of what happens in position space as the temperature is varied. 



2. Predicted Signal in Heavy Ion Collisions 

Let us consider the flux tube scenario as the fireball of QGP in heavy ion experiments expands 
and cools. The cross-over region does not possess a single percolating flux tube network, but it 
can still contain many finite clusters of flux tubes. Since the flux is directed, an obvious feature of 
every cluster is that any neighbour of a vertex is an antivertex and vice versa (see Fig]3|(d)). This 
is a topological feature, rather immune to the details of QCD dynamics. Conservation of baryon 
number implies that the vertices can only be locally pair-produced or pair-annihilated. As the QGP 
hadronises, the flux tube clusters start breaking up. After the chemical freeze-out stage, there is no 
more production or annihilation of vertices; every vertex ends up in a baryon and every antivertex 
ends up in an antibaryon. In the absence of subsequent large-scale diffusion, the radial propagation 
of (anti)baryons preserves the geometric pattern of (anti)vertices present at the chemical freeze-out 
stage, and the angular positions of the (anti)baryons seen in the detector can be backtracked to the 
angular positions of the (anti)vertices at the chemical freeze-out stage. This pattern of vertices can 
then be analysed for correlations and fluctuations, using techniques similar to those used to analyse 
the temperature correlations and fluctuations in the CMBR [^. 

Numerical simulations can calculate the equilibrium congelations between flux tube vertices 
in the 3-dim position space, which can then be projected onto the observed surface of the fireball. 
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Flux tube clusters are bipartite graphs of degree 3. Their alternating neighbour pattern of vertices 
and antivertices can yield two-point baryon number coiTclations similar to the oscillatory two-point 
charge congelations in ionic liquids, illustrated in Fig|I](b). The key common ingredient in the two 
cases is the hai^d-core repulsion between the objects involved. 

An oscillatory behaviour of the two-point baryon number coiTclations requires the following 
(in a mimicry of the Sakharov conditions): (1) A mechanism for local baryon-antibaryon pair pro- 
duction must exist. This is inherent in the QCD dynamics. (2) The system must have a non-zero 
chemical potential (analogous to a Fermi surface in condensed matter systems). This avoids spec- 



tral positivity constraints, as in case of symmetric systems [11], and is true in experiments. 
(3) The dynamical evolution must not be in equilibrium. This favours fragmentation of flux tube 
clusters during hadronisation while suppressing vertex-antivertex annihilation. Production of a 
sizeable number of antibaryons in experiments from an initial state that has none, confirms it. 

For concreteness, consider a homogeneous and isotropic fluid consisting of discrete objects 
with charges ^, = ±1 at positions r,. Then the average density Pq is independent of the position 
and the average two-point correlation function gq{fi,fj) depends only on r = |r; — fj\: 

Pq = {Y.^i^{ri)^ ^ pqgq{r) = {Y^qi5{r-\fi-rQ\)Y (2-1) 
i i/O 

Interactions fade away at long distances, and so g{r — )• oo) = 1. For objects with hard-core repulsion, 
g{Q) = 0, and beyond the hai^d core g{r) tends to its asymptotic value exhibiting damped oscillations 
JT2[]. In particular, the distance scale of the oscillations is determined by the inter-object separation, 
and the amplitude of the oscillations is determined by how tightly the objects are packed together. 

To quantify the correlation between vertices and antivertices, it is useful to compai^e the conge- 
lation functions g,, and ^|,,| on the same data sets, respectively obtained with and without the factors 
of qi in Eq.(^). They are sketched in Fig.|l|(b), and their comparison reduces systematic effects 
arising from varying fireball sizes. g\y\ (r) is rather insensitive to interactions other than the hard 
core. Its first peak is the most informative one regarding the equilibrium fluid properties — with the 
location, the height and the area under it respectively giving estimates of the inter-object spacing, 
the fluid compressibility and the number of nearest neighbours. On the contrary, successive cone- 
lated neighbours contribute with opposite signs to gv{r). For a system with no conelations between 
vertices, the probability of occunence of a vertex or an antivertex at any location is proportional to 
its overall density, and gv{r) behaves the same way as g\^,\ (r). Thus the contrast between gv{r) and 
g\v\ ix)^ which is maximum at the first peak, is a measure of the vertex-antivertex conelations. 

Note that a break-up of the flux tube clusters smoothens the oscillations of the two-point 
conelation functions, similar to what happens when a fluid is warmed. Still, if large enough flux 
tube clusters remain in the QGP, they would contribute to the contrast between gv{r) and (r). So 
the strength of the contrast measures the extent to which the nearby deconfinement phase transition 
influences QGP properties in the cross-over region. Specifically, the features to be quantified are 
the distance scale and the amplitude of the oscillatory signal. The former is expected to be the 
inter-baryon separation (~ 2 fm), while the latter indicates how tightly or softly the QGP is packed. 

Projection of the 3 -dim conelation functions onto the observed surface of the fireball smeai^s 
their oscillatory structure. That is an easily calculable effect as long as the inter-vertex separation 
is smaller than the radius of the fireball [^. More importantly, the angular coverage of the ex- 
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perimentally observed hadron distributions is limited. A partial Fourier expansion, based on the 
rotation symmetry around the beam axis and parity, helps in the search for congelations by orthog- 
onal separation of scales. Labeling the angular distributions by the unit vector h{d, 0), we have 

-| OO 1 CO 

K«) = ^I £ {b{h)b{n')) = ^Y, £ (2-2) 

y m (y=±m=-£>o (T=±m=-oo 

These expansions of baryon number distributions can be inverted as 

C^(0>0') = ^y^ d(^l d4' e'^^^-^\b^{h)b^{h')) , b^{h) = -{b{h)±b{-h)) . (2.3) 

The real and symmetric coefficient functions C^{d,d') contain all the information about the two- 
point correlation functions, and their experimental values can be tested against model predictions. 

The biggest gap between the theoretical formalism and the experimental data is the fact that 
the detectors observe protons and anti-protons but miss neutrons and anti-neutrons. One has to rely 
on the isospin symmetry to assume that the observed subset of (anti)protons provides a faithful 
characterisation of the total baryon number distribution. Corrections also need to be estimated for 
only approximate equilibration of the fireball, non-uniformity of the QGP caused by the elliptic 
flow, baryon number diffusion subsequent to hadronisation and development of the hard baryon 
core during hadronisation. Nevertheless, it is imperative to extract the two-point baryon number 
coiTclations from the experimental data and compare them to theoretical estimates. That would tell 
us a lot about where exactly the QGP lies between the extremes of a rigid crystal and a dilute gas. 
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